Bessel beam propagation: Energy localization and velocity 
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The propagation of a Bessel beam (or Bessel-X wave) is analyzed on the basis of a vectorial 
treatment. The electric and magnetic fields are obtained by considering a realistic situation able to 
generate that kind of scalar field. Specifically, we analyze the field due to a ring-shaped aperture 
over a metallic screen on which a linearly polarized plane wave impinges. On this basis, and in the 
far field approximation, we can obtain information about the propagation of energy flux and the 
velocity of the energy. 

The motion of a Bessel beam is of great interest in physics both for its characteristic as a non-diffracting beam 
US Hll SHU, and for its implications with regard to the topic of superluminality @ H E3, El EHH . Extended 
studies have been devoted to these subjects from both an experimental and a theoretical point of view. However, in 
spite of the many efforts devoted to this topic, no definite answer has been found about the amount of the energy 
transfer and its velocity. 

As far as Bessel beam propagation is concerned, the problem is mainly related to the difficulty in finding the 
vectorial field that describes this system, while, on the contrary, the field in the scalar approximation is well-known. 

The purpose of the present work is to investigate the propagation of a Bessel beam on the basis of a vectorial 
treatment. This kind of approach allows us to obtain information regarding the propagation of the energy flux and 
the energy mean velocity. 

A Bessel beam consists of a set of plane waves with directions of propagation s = a\i+ /3ij + 71k which makes 
the same angle 8q (0 < 9q < tt/2) with the z-axis (hence 71 = cos So for all the plane waves). Tn spherical coordinates 
(p, 9, p), the direction of propagation is specified by 

u\ — sin 0o cos if , (3\ = sin 9o sin ip 7i=cos#o- (1) 
Thus, for propagation in vacuum or air, each one of these waves, at the point x, y, z, can be written as 

u(P) = uod(p enp[iko(aix + ft\y + jiz)) exp(— iuit), (2) 

where Ugdp is the amplitude of the elementary wave, u> is the angular frequency, k$ — uj/c is the wavenumber, and 
x, y and z denotes the Cartesian coordinates of P. In cylindrical coordinates p, tp, z around the z-axis 

x = pcosip, y = psintp, z = z, 
and the total field U , given by the superposition of all the waves J2Jl, can be obtained by integrating over dip, that is, 

U = uq / ex.p[iko(aix + [3\y + jiz)] exp(— iuut) dp 
Jo 

= 2tt uq Jq (k p sin 9 q) exp (iuj— cos0 o ^ exp(—iujt) (3) 

where Jo denotes the zero-order Bessel function of first kind[l4|. 

The scalar field of Eq. (0 is known as a Bessel beam (or Bessel-X wave), the unusual features of which are 
- that it does not change its shape during propagation, since its amplitude is independent of z 15]; 
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FIG. 1: Scheme of the system considered for the theoretical analysis presented here. The metallic screen (over which the 
ring-shaped aperture is placed) and the converging system must be considered as being of infinite dimension. Details are shown 
in the inset. 



- that it propagates in the z direction with phase and group velocities v — c/ cosOq larger than c|3.l9L ITol[Tf^ |. 
Both the above mentioned characteristics can be analyzed in detail by means of a vectorial treatment, since the scalar 
field J3J represents an approximation of an electromagnetic field, and only a knowledge of the vectorial field (and of 
the Poynting vector in particular) can provide detailed information about the energy propagation. 

Vectorial fields with amplitude proportional to the zero-order Bessel function can be found in different ways. 
However, in order to derive just the vectorial field describing a system which has Eq. Q as scalar approximation, 
we have to consider a realistic situation that is able to generate a field of that kind. For this purpose, let us consider 
the system of Fig. ^ which consists of a ring-shaped aperture, of radius r, over a metallic screen on which a linearly 
polarized plane wave impinges (from the left). The ring is placed on the focal plane of the converging system C with 
focal length / 3> A, A being the wavelength. Let us consider the impinging electric filed to be polarized in the i 
direction, and the thickness d of the ring to be very small with respect to A. We may assume the element dA of the 
ring, at S{r, tp, —f), to behave as an elementary dipole, parallel to i, with amplitude proportional to dip (this requires 
a suitable choice of the transparency of the ring at S). Thus, the associated field at the optical center O of C has the 
well-known characteristics of the far field radiated by an elementary dipole, that is to be a spherical wave centered 
at S, with the electric field e in the meridional plane of the dipole through O, and perpendicular to the direction R 
from S to O. We can write 

R = — r cos ip i — r sin tp j + / k, 

e = e x i + e y j + e z k. (4) 
The two characteristics mentioned above can be written as 



e • R = — r cos ipe x — r sin ip e y + f e z = 

(R x e) • i = / e y + r sin tp e z = 0, (5) 

where we disregarded a phase factor exp(i/ci?). By solving Eqs. Q we obtain 

/ r sirup \ 

e y = -r cos ipe x — , , — (6) 



r 2 sin tp + f 

e z =rcosipe x [ . / — - J . (7) 
\ r z sin tp + j * J 

At this point, we can assume that the field emerging from the optical converging system is a plane wave propagating 
in the direction of R, with amplitude proportional to the amplitude at C of the incident field, namely 



e_pr< 
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FIG. 2: Electric fields E x , E y and E z normalized to 2neo exp(i£z) (continuous lines) vs. p, as given by Eqs. I121 - H14j . for 
k = 2, ij) = 10° and #o = 10°. In E z phase factor e 17V < 2 is disregarded. The dashed line represents the normalized scalar field 
© for the same parameter values. 
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(e x i + e y j + e z k) e~xp[ik(ax + f3y + 72:)] 
1 
Z 



(ai + (3] + 7k) x (e x i + e y j + e z k) exp[ik(ax + (3y + 72)] 



(8) 
(9) 



where a = — costp sin 6q, (3 = — sin ip sin #0 , 7 = cosf?o are the director cosines, Z is the free-space impedance, and 
the temporal factor exp(— iu>t) is omitted. The total electric field E will be given by the superposition of all edip 
contribution arising from the dipoles, and results in 



E 



edip = J (e x i + e y ] + e z k) exp — p (^— ^ cos(<p — ip) + cos^oz | dip 
exp (ik cos 9qz) / (e^i + e y j + e^k) exp [— ikp sin # cos(<p — ip)] dip. 



With reference to Eqs. © and J2J, it is expedient to choose 



£0 
P 



(r 2 sin 2 ip + f) 



(10) 



(11) 



this condition can be experimentally obtained by a suitable choice of the transparency of the ring as a function of 
(p. Thus, by substituting Eqs. (JBJ), Q) and into Eq. (|10|) . and by recalling that r = /tan6*o, we finally obtain 
(calculations are rather cumbersome but of no difficulty) 



E„ — 



En 



2Tre e liz I Jq{t]p) + tan 2 6 



Mvp) ~ 



2tt e e^ z 



sin 2ip 



tan 2 6» Mvp) 



VP 
VP 



cos z 4> Jo(vp) 



2Ji(vp) 
VP 



E z — — 27ri eoe'^ z [tan^o cos^Ji(?7p)] , 



(12) 

(13) 
(14) 



where £ = fccos#oi ?7 — fcsin# j an d J\ denotes the first-order Bessel function of first kind. 

Equation (|12|) (the main contribution of the electric field) describes a field different from the scalar field of Eq. © 
because of the presence of the term depending on tan^o. However, for 6q <C tt/2 (r <C /), as in the present case, this 
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term is negligible. We note that also the dependence on if) (which is absent in the scalar approximation) is negligible, 
and may be due to the approximation indicated in |l7j . In Fig. [2 we report the normalized value of E x , E y and E z vs 
p, for #o = 10°, together with the scalar field of Eq. The scalar field is practically coincident with E x . Therefore, 
we can conclude that the vectorial field derived above has Eq. @ as its scalar approximation, at least for r <C /. 
The magnetic field can be derived by Eq. (JSJ, and results in 



H x = (15) 
H y = ^eoe^-^-Mvp) (16) 

Zj COS Uq 

■ 2tt : P7 . , sin6> T , . 
Z cos z &o 

From a knowledge of the electric and magnetic fields, we are now in a position to evaluate the mean density of the 
energy flux which is defined as one half of the real part of the complex Poynting vector [Tsl Il9j 
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FIG. 3: Energy flux S z as given by Eq. 1)190 normalized to its maximum value vs p, for a few values of 6o- For the sake of 
completeness, in the inset S z is shown for larger values of p, even if only the region around the main peak is of physical interest. 
Parameter values are as in Fig. 



S = -Re(ExH*). (18) 
For the fields lfT21) - (|Tl|l and ljT5|) - (|T7Jl . it turns out that S has only the fc-component 

S 2 = \ (E X H*) , (19) 

that is the propagation of the energy flux occurs only in the z-direction, in accordance with the information given by 
the scalar field J3J. Moreover, since the flux is independent of z, the energy propagates with no deformation. 

In Fig. |21 the behavior of the energy flux (|19fl is shown as a function of p for a few values of #o- We note that, 
for #o very small (nearly plane wave) the flux is nearly independent of p, while when the beam originates the flux 
increases by increasing and tends to concentrate near p = 0, that is, along the z-axis. Thus, for small values of p 
(that is, in the proximity of the z-axis), the power supplied by a Bessel beam is always greater than the one due to a 
plane wave. 

As for the velocity of the energy, from Eq. (0J it follows that in the scalar approximation the dependence of the 
field on t and z occurs only through the quantity (z/c) cos#o — t and, therefore, the field propagates with velocity 
v = cj cos 6*o. On the basis of these arguments, it could be concluded that also the energy propagates with a velocity 
v e greater than c. In the vectorial treatment we can evaluate the energy velocity v e as |l8l Il9| 
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FIG. 4: Energy velocity normalized to the light velocity vs p, for k = 2, ip = 10°, So = 10° (blue line) and for 8o = 12° (red 
line). Dashed lines represent the energy flux relative to the same parameter values. 



\{eE-E*+pH-H*)' V 1 

where the quantity (1/4) (eE ■ E* + pH ■ H*) is the total mean density of energy which can be evaluated with the help 
of Eqs. (0-031) and (fT%|) - fT7|) . 

In Fig. 21 w e report the normalized velocity of energy as a function of the radial coordinate p, for 8q = 10°. The 
velocity is found to be equal to c from p = up to near to the first zero of the Bessel function: that is, the beam 
moves like an almost rigid system, in spite of its dependence on p and In the proximity of the first zero of the 
Bessel function, the velocity decreases and tends to zero. Naturally, the zero in the velocity does not represent a stop 
of the motion but, more simply, the absence of energy flux. In this situation, the concept of velocity has no physical 
meaning. 

Some remarks must be made on this surprising result. In fact, we recall that for propagation in vacuum "if an 

energy density is associated with the magnitude of the wave the transport of energy occurs with the group 

velocity, since that is the rate of which the pulse travel along" 20J. If the definition of the energy velocity as given 
by Eq. H20[) is applicable also to a Bessel beam (or, more generally, to localized waves), it is not clear what kind of 
physical mechanism makes the energy velocity different from the phase and group ones. 

We wish to recall that the present analysis was performed for an ideal system in the far field approximation. For 
a real system, we have to take into account the finite dimension of the converging system, which limits the field 
depth and introduces diffractive effects. The role of diffraction, together with the analysis of the near field (as in real 
experimental situations), make the problem much more complicated, and is beyond the purpose of this paper. 
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